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Abstract 

Recent work has raised considerable interest on the nature of thin metalhc 
wires. We have investigated the melting behavior of thin cylindrical Pb wires 
with the axis along a (110) direction, using molecular dynamics and a well- 
tested many-body potential. We find that — in analogy with cluster melting — 
the melting temperature Tm{R) of a wire with radius R is lower than that of a 
bulk solid, T^, by Tm{R) = Tj^ — c/R. Surface melting effects, with formation 
of a thin skin of highly diffusive atoms at the wire surface, is observed. The 
diffusivity is lower where the wire surface has a flat, local (111) orientation, 
and higher at (110) and (100) rounded areas. The possible relevance to recent 
results on non-rupturizing thin necks between an STM tip and a warm surface 
is addressed. 

PACS numbers:64.70.Dv,82.65.Dp,61.46.+w 
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It is well known — experimentally and theoretically — that small atomic clusters have a 
melting temperature which is lower than that, T^, of the infinite bulk solid of the same 
element The main physical reason for this is that the liquid- vapor interface free energy 
7to is generally lower than the average solid-vapor interface free energy jsv Therefore, when 
the particle size is decreased and its surface to volume ratio increases, the "bulk" liquid free 
energy at Tm must become higher than the "bulk" solid free energy to compensate for the 
advantage of the liquid at the surface, and Tm must decrease. 

Thin wires with a microscopic diameter should exhibit a similar melting behavior — 
although with significant differences due to their cylindrical symmetry and infinite extension 
along the axis — but they have not been extensively investigated so far, probably due to lack 
of practical motivations. However, considerable interest in the thermal behavior of thin wires 
has been recently raised by some Scanning Tunneling Microscopy (STM) experiments. 
Abrupt changes of the tunnel resistance have been observed when the tip-surface separation 
distance becomes sufficiently small, indicating jump into contact [Q. Molecular dynamics 
(MD) computer simulation studies [^-0 have indicated that contact is often realized by 
means of a narrow "neck" connecting the tip with the sample surface. In an MD study of 
a Au tip interacting with a Pb(llO) surface at high temperature this neck is seen to be 
liquid in nature in its thinner portion, even when the temperature is considerably lower than 
Tm- Moreover, the jump to contact between tip and surface via such a partly liquid neck 
was predicted, because of this, to occur more readily as temperature was raised. Recent 
STM experiments by Kuipers et al. also performed on Pb(llO) using a W tip over a 

temperature range from about 0.5 T/y/ to surface melting, indicate that the neck elongates 
up to considerable lengths (several thousands of A) , depending on the retraction speed and 
on temperature, and then it breaks. To shed some light on the behavior of such a thin and 
long neck, we have started an investigation of cylindrical metallic wires. 

(i) Simple theory: For spherical particles of radius i?, a melting temperature Tm{R) 
is obtained phenomenologically by equating the Gibbs free energies of solid and liquid 
spherical clusters, assuming constant pressure conditions: 
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where ps and pi are the sohd and hquid densities, is the bulk latent heat of melting and 
•jsv and are the solid-vapor and liquid-vapor interface energies respectively. The surface 
energy anisotropy of the solid is not taken into account, and the possibility of inhomogeneous 
phases (such as a liquid skin wetting the solid cluster) is also neglected. In spite of these 
approximations, the 1/R behavior is approximately correct for clusters of sufficiently large 
size, although in some systems Tm drops faster than predicted ||l|, [lO|Jl3| . In the case of Pb, 



which is the subject of our MD study, deviations from ~ 1/R appear to be small even 
at small sizes |1TT|JT^ , For a wire, following a similar procedure, we equate the Gibbs free 
energies per unit length of the solid and liquid wires at constant pressure and temperature, 
Gg = Nps + 27ri?7st, and Gi = Npi + 27ri?7i„, where N is the number of atoms per unit 
length, Ps and pi are the chemical potentials of the solid and liquid phases respectively, and 
R is the wire radius. Since near T^, N{ps — pi) = VpsL^(T — T^)/T^, we obtain for the 
melting temperature Tm{R) of a thin wire 

Tl - T^R) , ^ / (2) 
n, p.L''R\'"' \p,) ^•") 

Hence, the melting temperature of a wire should be depressed, approximately, by half the 
corresponding depression of a spherical cluster. In this simple model the latent heat of 
melting per atom (averaged on all atoms) decreases with exactly the same law when the size 
is decreased: 

These relations imply the existence of a critical radius Rc = {'jsv — i^Y^'^liv) / (PsL^) corre- 
sponding to Tm{Rc) = 0, L{Rc) = 0, below which the solid wire is not stable. 

(ii) MD Simulation: For our simulation we have chosen lead, and a (110) orientation 
for the wires. Pb was a natural choice as the melting properties of its surfaces and clus- 
ters have been extensively investigated both experimentally and theoretically. Interestingly, 



Pb(lll), Pb(lOO) and Pb(llO) exhibit different melting properties, that is, respectively 



non- melting, incomplete melting and surface melting |T^-|T6[. On small Pb particles, the 
formation of a liquid skin as a premelting effect has been observed by microscopy experi- 
ments |Tl|]. We have used a "glue" many-body potential, already tested in previous surface 



JT7|JT^ and cluster studies. In particular, the different melting behavior of the main 



surfaces of Pb was found, in excellent agreement with experimental results. Because the 
model potential is untested, and therefore unreliable at really low coordination numbers, 
we restrict our study to i? > 10 A. By choosing a (110) wire axis, (111), (100) and (110) 
orientations are simultaneously present on the cylindrical surface of the wire. The different 
melting behavior of these facets can make the details of premelting phenomena particularly 
interesting. Moreover, Pb(llO) is the system on which the STM experiments discussed above 
have been done, and one could conjecture that the neck might retain the same (110) 
orientation of the sample surface. 

We have performed several MD runs at different temperatures and for different wire 
radii. All samples are made up by 36 (110) layers, repeated along the wire axis z by periodic 
boundary conditions. The initial configurations are prepared by including all atoms in an 
fee lattice whose distance from the axis is smaller than an assigned radius. At T = i^, 
well defined (100), (110), and (111) facets are present (Fig. 1). This roughly cylindrical 
geometry is expected to be relatively close to the equilibrium situation, since the surface 
energy anisotropy of Pb is known to be rather small [0 . For the system sizes investigated, 
this procedure typically results in octagonal cross sections. We have studied nine samples, 
with a total number of atoms = 1134, 1206, 1314, 1710, 1854, 2286, 2574, 3294 and 4230, 
with radii between i? = 13 A and i? = 25 A. 

All MD runs were performed at constant temperature and constant wire length, where 
the box length along z is adjusted according to the bulk thermal expansion curve. While 
it would also seem possible to operate at constant axial stress and varying length, it is not 
feasible to achieve stability under such conditions when the wire becomes liquid. The time 
step is At = 7.32 x 10^^^ sec. The initial samples with atoms in ideal fee positions, were 



first relaxed by simple quenching. Then, an annealing cycle at room temperature, followed 
by a quench to OK, is done to check the stability of the original structure. In most cases, due 
to creation of defects during the annealing, the final energy was slightly higher than that of 
the previous relaxed structure. The absence of any major structural changes indicates that 
the original structure was optimal or nearly optimal. Each system was then heated by MD 
runs of 20000 time steps (147 ps) with a temperature step of 50K up to 300K and of 25K 
above 300K. In proximity of Tm we further reduced the temperature step to 5K and used 
longer equilibration runs of 40000 time steps (293 ps) at each temperature. We monitored 
several quantities, such as the internal energy E, and the diffusivity as a function of T. 

The caloric curves E{T) for all the wires investigated are shown in Fig. |^. The dashed 
vertical line indicates the bulk melting temperature of our model determined from bulk MD 
runs at zero pressure, =618K (to be compared with 600. 7K of real Pb). In each curve 
one can identify four regions, labeled as A, B, C, D in Fig. |^. In region A the wire is 
solid, and the slope corresponds to the Dulong- Petit specific heat. Visual examination of 
the samples confirms that diffusion is non-existent, or very low in this region. In region B, 
E{T) exhibits an upward curvature, where the specific heat gradually increases by a small 
fraction of the total heat of melting. This increase is associated with loss of solid rigidity, 
and onset of diffusion at the wire surface. For example, the wire cross-section in Fig. |l|c, 
taken at T = 570K for a wire with R = 22.5 A, clearly shows that melting starts first at the 
wire surface, while the inner region is still ordered. The surface diffusivity is anisotropic and 
appears to be lower in correspondence with (111) facets. This is reasonable in view of the 
known non-melting behaviour of Pb(lll) ||T6| , p!7| . The lateral view (Fig. ||d) also shows the 
presence of (111) crystalline patches on the wire surface. Presumably, the wire size is too 
small to exhibit complete crystallinity of (111) facets, which should be expected in the limit 
of macroscopic sizes. (100) patches are also sometimes observed. (110) faces appear to be 
the first to melt with increasing T, followed by (100) and then (111). Another interesting 
feature showed by Fig. |I|c is that (110) facets — and, to a lesser extent, (100) faces — have 
shrunken in size with respect to the low-temperature shape, and are considerably rounded. 
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A rounded profile is characteristic of a rougliened surface, where steps can proliferate freely. 
This observation is therefore qualitatively consistent with the experimental observation of a 
roughening transition of Pb(llO) around 415K |[T^, and is further corroborated by surface 



simulations results |]I8[. Region C consists of an abrupt nearly vertical jump corresponding 
to the melting transition of the whole wire, and its temperature and magnitude are discussed 
below. In region D the wire is completely liquid, and E{T) grows again linearly. 

It is interesting to note that the internal energy jump (region C), corresponding to the 
latent heat, remains sharp and discontinuous: the wire suddenly melts when its melting 
temperature is reached, regardless of the quantity of solid at the surface which is already 
melted. This sharp jump yields a rather precise value of for all the sizes investigated. 
There is a clear decrease of with decreasing wire radius (Fig. To verify the predicted 
approximately factor of two of between Eqs. (|l]) and (0) we compare the behavior of wires 
with that of Pb clusters investigated by Lim et al. [0 using the same force model. Fig- 
ure ^reports as a function oil/R for wires (solid circles) and clusters (open squares). 
The dashed lines correspond to the simple thermodynamical models of Eqs. (|T]) and (^, 
using = 544mJ/m2, = 476mJ/m2, ps = 0.0321 atom/A^, pi = 0.031 atom/A^ and 

= 0.0494 eV/atom. For both geometries, the 1/R behavior is followed with fairly good 
accuracy, as well as the approximate ratio 2 between the slopes. 

The energy jump on melting corresponds to the latent heat L{R). As a result of the 
premelting effects discussed above, this jump is slightly smaller than the latent heat that 
would be obtained if the system were entirely crystalline just below T^. In order to verify 
the prediction of Eq. (^, where premelting effect are not considered, we have constructed a 
modified latent heat L*{R), defined as the difference between the energy of the liquid and 
the energy of a hypothetical complete solid at T^, obtained by ignoring surface melting, 
and extrapolating linearly E{T) from its low-temperature behavior. L*{R) normalized by 
bulk latent heat, Lf, = 0.056 eV/atom determined from our model, is plotted in Fig. ^ as a 
function of l/i? (solid circles), together with the data relative to spherical clusters from ref. 
jT^ (open squares) and with Eq. (|^) (dashed line). For both geometries the 1/R behavior 
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is again approximately true, although not quite as accurate as for T^. However, deviations 
from Eq. are more noticeable for clusters than for wires, and for all the smaller sizes. 
A faster decrease of L compared with Tm with decreasing size has been also observed in 
previous Au clusters simulations |jl3[. These deviations can be tentatively attributed to the 
reduced entropy of atoms at the liquid metal surface compared with those in the bulk liquid. 

Our results on wire melting may be of relevance to the neck between an STM tip and 
a metal surface. While a thick crystalline neck can still be slowly pulled out of a surface 
via tip-induced diffusion, as explained in Refs. and it is possible, as indicated by 
MD simulations that at least a part of the neck could be extremely thin, R ~ SA. Our 
extrapolated melting temperature for a Pb wire so very thin, is lower than room temperature. 
Hence, that neck portion would be completely liquid above 300K, allowing matter to flow 
freely along its axis even under fast traction. The possible existance of these ultra-thin 
necks might be checked by looking for, e.g. transversally quantized electronic levels and 
related jumps in the IV characteristics, and their liquid nature could be of importance for 
friction problems. 

In conclusion, thin wires, like small clusters, melt at a temperature Tm{R) which is lower 
than the bulk melting temperature, and is found to be approximately a decreasing linear 
function of inverse radius. The latent heat follows a similar law, but drops slightly faster 
with decreasing size. Wire melting is preceded by surface melting effects of its outer skin, 
similar to those exhibited by fiat surfaces and by spherical clusters. Possible relevance to 
the nature of thin tip-surface necks has been pointed out. 
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FIGURES 

FIG. 1. Snapshot views of the MD sample with N = 3294 and R ~ 22.5 A. a) and b) is 
the initial sample at OK. a) is the top view of the wire cross section. Wire axis (z-direction) is 
perpendicular to this cross section, b) is the side view from an angle such that all different laterals 
are seen, c) and d) is corresponding views at a temperature just below the melting temperature. 
The grayness of each atom is proportional to its square displacement during a MD run. White 
indicates a non-diffusing, crystalline behavior. 

FIG. 2. Internal energy as a function of temperature is shown for different wire samples. Num- 
bers labeling the curves are the number of atoms in MD cell. Dashed vertical line indicates the 

bulk melting temperature for lead determined from MD. Discontinuity in the curves shows the 
melting transition. 

FIG. 3. Melting temperature of different clusters as a function of inverse radius. Solid circles 
refer to wires, open squares to spherical clusters, the dashed lines are the predictions of the simple 
phenomenological theory (see text). 

FIG. 4. Latent heat of melting for different clusters as a function of inverse radius. Solid circles 
refer to wires, open squares to spherical clusters, the dashed lines are the predictions of the simple 
phenomenological theory (see text). 
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